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Abstract 



Consider a d-diniensional antiferromagnet with a quantum disordered ground 
0^ ■ state and a gap to bosonic excitations with non-zero spin. In a finite exter- 

nal magnetic field, this antiferromagnet will undergo a phase transition to a 
ground state with non-zero magnetization, describable as the condensation of 
Q . a dilute gas of bosons. The finite temperature properties of the Bose gas in 

^" \ the vicinity of this transition are argued to obey a hypothesis of zero scale- 

factor universality for d < 2, with logarithmic violations in d = 2. Scaling 



f^ ■ properties of various experimental observables are computed in an expansion 

^J . in e = 2 — d, and exactly in d = 1. 
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I. INTRODUCTION 

Numerous experiments |I| have by now examined the properties of 5 = 1 spin chain 
antiferromagnets which posses the Haldane gap 0. More recently, these antiferromagnets 
have been placed in a strong magnetic field 0, and found to display evidence of a zero 
temperature phase transition to a state with a non-zero ground state magnetization. Theo- 
retical studies 1^-0 have also examined this transition at zero temperature, and developed 
a picture of it as the Bose condensation of magnons with azimuthal spin 5"^ = 1. Although 
this possibility has not been considered before, it is not difficult to see that the condensation 
of magnons in a finite field should occur in quantum-disordered antiferromagnets in any di- 
mension d (provided, of course, the magnons continue to behave as bosons). In this paper 
we shall present a general theory, in dimensions d < 2, of the finite temperature properties 
of quantum antiferromagnets in the vicinity of such a field-induced quantum transition. 

We begin by elucidating the precise conditions under which our results apply. Consider 
an antiferromagnet with a quantum disordered ground state. The Hamiltonian must posses 
at least an axial symmetry i.e. at least one component (say z) of the total spin must 
commute with the Hamiltonian. The lowest excitation with non-zero spin must be separated 
from the ground state by a gap, and behave like a quasiparticle with bosonic statistics. 
These conditions are clearly satisfied by Haldane gap antiferromagnets in which the in- 
plane anisotropy can be neglected (the compound NENP does have a rather small in-plane 
anisotropy [Q], and the restrictions this places on applying our results to on experiments in 
NENP will be discussed later). In d = 2, the 5* = 1/2 kagome Heisenberg antiferromagnet 
is probably the most accessible candidate upon which our results can be tested - it has been 
argued that the magnons in this system are bosons |^,0. 

Now place this antiferromagnet in a magnetic field pointing along the direction of axial 
symmetry. The eigenenergy, e{k) of single magnon (boson) quasiparticles with momentum 
A; in a field H then takes the form: 

e{k)=/\ + ^-giiBS.H (1.1) 

Here A is the magnon gap and m is the quasiparticle mass, both determined in the zero field 
antiferromagnet. For the S = 1 Haldane gap chain, we have the azimuthal spin Sz = 1. In 
the kagome antiferromagnet, the magnon excitations have been argued to be spinous [§0: 
we therefore expect Sz = 1/2. Later, we will discuss how the field induced transition in 
the kagome antiferromagnet may offer a way of experimentally/numerically determining the 
value of Sz- 

The antiferromagnet will undergo a T = field-induced transition at the field H = H,, 
which is given exactly by 

gi^BSzH, = A (1.2) 

In the vicinity of this transition, we may describe low energy properties of the antiferromag- 
net by just studying an effective Hamiltonian for the bosonic magnons IQ. The remainder 
of this paper will therefore consider properties of the following coherent state path integral 
over the Bose field \E'(x,r) 
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where x is the d-component spatial co-ordinate, r is the Matsubara time, the chemical 
potential 



/i = g^isSzH - A, 
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This field theory has a phase 
13 • They 
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and f is a repulsive interaction of a short range ~ A" 

transition exactly at T = 0, /z = which has been studied by Fisher et. 

identified the upper-critical dimension as c/ = 2, above which the interaction v is irrelevant. 

For c? < 2, f is relevant, but the exponents were nevertheless found to have trivial values: the 

dynamic exponent z = 2, the correlation length exponent v = 1/2, and the field anomalous 

dimension r] = 0. The triviality of the exponents is partially related to the fact that the 

parameter tuning the system through the transition, fi, couples to a conserved quantity - 

the density of bosons; any such transition [|l^] must have zu = 1. The structure of the 



d < 2, finite-f fixed point is thus very unusual: despite describing a non-trivial, interacting, 
critical field theory, the exponents associated with all the relevant directions away from 
this fixed point are trivial. In this paper, we shall show that the interactions are crucial 
in determining the finite temperature properties of the Bose gas near this fixed point. The 
fixed-point interactions are needed to preserve hyperscaling for d < 2 and lead to highly 
non-trivial scaling functions for the finite-temperature correlations. 

Before we state our zero scale-factor universality hypothesis for Z in its most general 
form, it is helpful to consider one of its simple consequences at T = 0. Examine the ground 
state boson density n = (|\l/(x, t)P) as a function of /x for small /x. This problem was studied 
many years ago for the d = 3 hard-sphere Bose gas ||lT| with the result 
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where a is the hard-sphere radius, and 6{x) is the unit step function. Note that, in addition 
to its dependence on m and /i, the boson density is sensitive to the nature of the boson-boson 
interactions (measured by the hard-sphere radius a). A different choice for the boson-boson 
repulsion would lead to different result for n. The situation in dimensions c? < 2 is however 
strikingly different; one manifestation of the zero scale-factor universality is that for small fi 



n-i^fcOi,) 



rf < 2, T = 
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where C is a universal number i.e. independent of the details of the interactions between the 
bosons. We will determine Cma.d = 2 — e expansion; its exact value in ci = 1 is known M : 



C = l/vr. The universality of C is a direct consequence of having a finite-couphng fixed 
point describing the onset at /z = 0. For d > 2, interactions are irrelevant, which leads to a 
violation of hyperscaling, and a dependence of the density on the nature of the microscopic 
interactions (as in the a dependence of ( |1.5|) ). Further, n will depend linearly on /i for all 
d > 2. Precisely in d = 2, as we shall see below, (|1.6|) is violated only by a logarithmic 
dependence on the microscopic interactions. 

For d < 2, the combination of the presence of hyperscaling, and the absence of any 
anomalous exponents in the leading critical behavior, leads to remarkably universal finite 
temperature properties. As in Ref. jT^, we may use finite-size scaling to deduce scaling 
forms away from /i = 0, T = 0. However the absence of any anomalous dimensions [z = 2, 
rj = 0, u = 1/2) means that the usual two scale-factor universality JTSf is now reduced to a 



zero scale-factor universality, described more precisely in the following subsection. 

A. Zero scale- factor universality 

In simple terms, this universality is just the statement that all response functions are 
universal functions of the bare coupling constants /i and m. There are no non-universal 



amplitudes; the usual case has two non- universal amplitudes, or scale- factors [|T3|. The 



universality can be stated more precisely in terms of the the boson Green's function 

G(x,r) = (TvI/(x,r)vI/*(0,0)) (1.7) 

where T is the ordering symbol in imaginary time r. After Fourier transformation as per 

f fh/(kBT) 

G{k,iun) = d^x cire-*(^-^-^"")G'(x, r) (1.8) 

this yields G{k,iujn) at the Matsubara frequencies Un = 27mT/h, from which the retarded 
Green's function G^{k, u) can be obtained by analytic continuation to real frequencies using 

G^{k, uj) = -G{k, iujn = uj). (1.9) 

The Lehman spectral representation of the Green's function implies that —ulm.G^{k, u) > 0, 
but does not constrain G^ to be an odd function of u. Our central result is the zero scale- 
factor universality of G^, which is equivalent to the scaling form 

where ^ is a highly non-trivial, but completely universal complex-valued function; naturally, 
A is independent of the nature of the boson-boson repulsion. An important property of A 
is that it is analytic at all finite, real, values of all three arguments. Similar scaling forms 
hold for other correlators of ^ - a particularly instructive observable is the local Green's 
function, Gf 

Gf{uj) = J-0^^G^{k,u) (1.11) 



If we take the imaginary part of this equation, it is expected that the resulting on-shell 
contributions on the right-hand-side will occur only at small momenta (determined by the 
frequency u), and the momentum integral is ultraviolet convergent. We may therefore 
deduce from (|1.1CI|) the scaling form 



ImGfiu) = -- 
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where F is a fully universal, dimensionless, positive function. It is quite instructive to 
consider the limiting behavior of F for small and large frequencies. We expect that Gf 
should be analytic at tu = at any finite T; this combined with the positivity condition on 
the spectral weight noted above, implies ImGf (a;) ~ uo for small uo at finite T. Therefore, 
from ( |1.12|) the scaling function, F, must satisfy F(ZJ,t) ~ |Up~'^/^ at small cJ (we use here 
and henceforth the notation uJ = hcu / (ksT) , and t = ix/{kBT)); the coefficient of the |cJp~'^/^ 
term is quite difficult to determine, and will be obtained in this paper only in a special limit. 
For large u, or short times, G^ should display essentially free particle behavior, as the dilute 
bosons have not had enough time to interact with each other. Using the free boson spectral 
weight we can deduce that 

F{^,t)= \(^^^2^ ^M as|a7|-.oo (1.13) 

Let us conclude this subsection by noting the precise conditions under which the system 
is in the critical region and (|1.10| ) and ( p..l2| ) are valid. We must have 



HkBT<t: ^^,viO) (1.14) 

Further the measurement wavevectors must satisfy 

fc<A (1.15) 

In d = 2 the zero-scale-factor universality is violated by a logarithmic dependence on 
the microscopic interactions. Furthermore, the scaling function A will have a singularity 
associated with the finite temperature Kosterlitz-Thouless transition. Zero scale-factor uni- 
versality does not hold for d > 2. 

B. Neutron scattering 

The dynamic information contained in G^ is directly observable in neutron scattering 
experiments. In Appendix |^ we discuss the relationship between the correlators of \1/ and 
antiferromagnetic correlations measured by the neutrons; this discussion is limited to the 
case where the quantum disordered phase has confined spinous i.e. S^ = 1, as in Haldane 
gap antiferromagnets. The relationship for the case of deconfined spinous {Sz = 1/2) is 
quite different WM and will not be considered in this paper explicitly. 



In the following, hcu will measure the energy lost by the neutrons in their interaction 
with the antiferromagnet. Consider first a scattering even in which the antiferromagnet 



undergoes a A5'^ = +1 transition. Then, from the discussion in Appendix ^ and the 
fluctuation-dissipation theorem we may conclude that (for the scattering cross-section is 
given by the dynamic structure factor S+_{k,uj) where 

-2ZlmG^ik,u;) 
S+-ik,uj)= ^_^_,^/,^^ (1.16) 

where Z is a non-universal quasi-particle renormalization factor between the magnon oper- 
ators and the ones that couple to the neutrons. (The wavevector k on the left-hand-side is 
measured from the antiferromagnetic ordering wavevector) The scaling result for S'+_ thus 
follows directly from ( |1.10| ). Next, consider scattering with ASz = —1 for the antiferro mag- 
net. The associated dynamic structure factor is then S-+{k,u) which is 

2ZlmG^ik, -u) 
S-+{k,uj)= ^_^.nu./k^T (1-17) 

By not resolving the energy of the scattered neutron, it is also possible to measure equal- 
time correlations embodied in the static structure factors. However the scaling properties 
of these observables are subtle, and require more careful interpretation. By performing a 
weighted frequency integral over the scaling limit of G^, we are implicitly only sensitive 
to frequencies much smaller than a high frequency cutoff like A/2. Thus in the following, 
our equal-time structure factors actually refer to scattering experiments in which energy 
transfers greater than A/2 are not integrated over. Not blocking out these events, will 
produce a background structure factor which may (as in the Haldane gap region defined 
below) overwhelm the universal part of the static structure factor which is considered here. 
Keeping this caveat in mind, we define the structure factor S-^—{k) 

S^-{k)= ^^S^.{k,u) (1.18) 

J A/2 2tT 

where the result is not sensitive to the precise locations of the limits. The correlator S |_ 

can be defined analogously. From (|1.18|) and ( p..lO|) we can deduce the scaling result 



S+.{k) = ZB+. { , ^^ , -^\ (1.19) 

where i3+_ is a universal scaling function determined completely by A (similarly for S ^{k)). 

Finally, using the fact that ^ and \1'* are canonically conjugate fields, it possible to deduce 
a frequency sum-rule on ImG^ which leads to 

S+_(r,t)=i3_+(r,t) + l (1.20) 



where we will henceforth use r = hk/ ^j2mkBT . 

We reiterate that all of the above results in this subsection refer only to the case of 
antiferromagnets with confined spinous. 



C. Uniform magnetization 

A second useful set of observables are those associated with magnetic fluctuations around 
k = 0. The uniform magnetization is a conserved quantity and, as a result, the quasi-particle 
renormalization factor Z does not appear in their scaling forms. The simplest of these is 
mean-value of the magnetization density M itself, which is of course related to magnon 
density by 

M = gfXBS.n. (1.21) 

This relationship, and the considerations of this subsection, apply to both confined {S^ = 1) 
and deconfined {S^ = 1/2) spinous. The mean-value of the magnetization is therefore 

M = gfiBS^G{x = 0,r = 0-) (1.22) 

and obeys the scaling form 

where the function JH is again dependent only on A. The result ( |1.6|) follows from ( |1.23|) . 



D. Phase diagram 

The above discussion shows that the scaling functions of a large number of experimental 
observables can be obtained directly from the primary scaling function A. The remainder 
of the paper is therefore devoted primarily to describing A and associated scaling functions 
in different parameter regimes. It is convenient to discuss the properties of A separately in 
three distinct regimes, which are analogous to those found by Chakravarty et. al. [|I5| in 
the d = 2 0(3) sigma model. These regimes are illustrated in the phase diagrams in Figs |l| 
{d = 1) and ^ (d = 2). The crossover boundaries between the regimes are delineated by the 
value of the dimensionless ratio fi/ksT (upto logarithmic terms in d = 2). We consider the 
three cases separately 



i. /i < -keT 



This is the analog of the quantum-disordered regime of Refs |T5|JT^ . Only a dilute gas of 
thermally-excited bosons is present, and their mutual interactions are weak. Properties of 
the quantum antiferromagnet can be described by a low-magnon-density expansion about 
the quantum disordered ground state. In d = 1, we identify this as the Haldane-gap regime 
in Fig |. 



2. |/i| < ksT 

This is the quantum-critical regime, in which z = 2 critical fluctuations are quenched 
in a universal way by the temperature. The value k^T is the dominant energy scale and 
universally determines everything: the boson density, spectrum, and interactions. The only 
small parameter which may be used to determine scaling functions is e = (i — 2. As we will 
see, this is not particularly effective in c? = 1 where, fortunately, exact methods are available. 

It is instructive to consider the physics of this region as a function of the measurement 
frequency u) (see Fig |). At large frequencies, uj ^ k^T ^ or short times, the effects flnite 
temperature have not yet become manifest, and the system displays the physics of the 
/i = T = critical fleld theory i.e. it is a dilute gas of bosonic quasiparticles with repulsive 
interactions. As one lowers the frequency through ksT /% there is a crossover to a novel z = 2, 
quantum-relaxational regime (Fig |^). Now each boson interacts strongly with thermally 
excited partners, leading to strong dissipation and overdamped quasiparticles. 

5. /i » ksT 

The behavior in c? = 1 and d = 2 is quite different and we will therefore consider the two 
cases separately. 

In d = 1, the ground state is a Luttinger liquid, which is itself a critical phase with 
z = 1 (Fig |I|). Again, consider the physics as a function of uj (Fig ^. For sufficiently 
large uj {fiuj ^ /i) we have the dilute bose gas physics of the ^ = T = critical point, 
similar to that discussed above for the quantum-critical region. At smaller u, there is a 
crossover (near hu = /i) to a Luttinger liquid-like region where we may as well assume that 
T = 0. However at small enough frequencies, huj ~ fc^T, the effects of a flnite temperature 
finally become apparent. The massless modes of the Luttinger liquid are then quenched 
into a, z = 1 quantum-relaxational regime, rather similar to the z = 2 quantum relaxational 
regime discussed above. This last crossover is, strictly speaking, a property of the z = 1 
critical point on the Luttinger liquid fixed line determined by the value of /i, and not a 
property of the z = 2 critical end-point at /i = T = 0. It is thus described by a reduced 
Luttinger liquid scaling function Al- We will obtain exact results for Al using an argument 
based on conformal invariance. We will also discuss an important compatibility condition 
between the scaling functions A and Al-, and show how the more general A collapses into 
the small /i limit of Al- 

In d = 2, (Fig|^) the ground state is a boson superfiuid, which survives at finite tempera- 
ture. There is a Kosterlitz-Thouless phase transition at a finite T, which has been studied in 



some detail earlier [0,^. In the superfiuid phase therefore, one has a large frequency dilute 
Bose gas behavior crossing over to a small frequency Goldstone phase with quasi-long-range 
order (Fig ^. 

The outline of the rest of the paper is as follows. We will begin, in Section ||, by a 
renormalization group analysis of the partition function Z. This will allow us to demonstrate 
the logarithmic corrections to scaling m d = 2 and obtain the leading terms of the in a e = 
2 — d expansion. These leading terms are consistent with the zero scale-factor universality. 
We will then, in Section |ITI|, turn to a discussion of the exact properties of A and Al in 
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d = 1. A brief summary and a discussion of relevance to experiments appear in Section |V. 
Three appendixes contain discussions of some peripheral points. 

II. RENORMALIZATION GROUP ANALYSIS 

The momentum shell renormalization group equations for Z have already been obtained 
by Fisher and Hohenberg ||r^. However their analysis of the equations was restricted to 
d = 2, fi > 0, and temperatures at or below the Kosterlitz Thouless transition. In this 
section we will extend this analysis to cover the remaining regions in d = 2 (Fig ^, and 
to dimensions (i<2inane = 2 — rf expansion (Fig p. The analysis [1^ proceeds by 
introducing an upper cutoff, A in momentum space, and replacing f by a contact interaction 
u = A~'^f (0). Degrees of freedom in a shell between A and Ae~^ are integrated out, followed 
by a rescaling of co-ordinates and field variables 

x' = e-'x t' = e-^'r (2.1) 

^'(x',r')=e'''/2^(x,r) 

Note that the there is no anomalous dimension in the rescaling factor for \E' as r^ = 0. 
Further, the scaling dimension of |\E'p is exactly d, as it must be for any conserved charge 
density [|l^ . 



It is convenient to consider the cases T = and T > separately: 

A. r = o 



The renormalization group equations are |]T^|TO 

^ = 2, (2.2) 

-=(2-d). i^u- (2.3) 

where K^ = Sd/{2-nY and Sd is the surface area of a unit sphere in d dimensions. We 
integrate these equations to a scale e^* where the system is non-critical i.e. when /i ~ 
ah'^h?' /{2m) where a is significantly smaller than unity, but not so small that the system is 
still critical. For /i > 0, the magnetization, M, and the boson density n = M/^gfisSz) are 
then given by 

^^^-.r/fgl (2.4) 

u{i*) ^ ' 

For (i < 2, and e = 2 — d small, u approaches its fixed point value 

« = ,^ w 9 g 2.5 



Inserting this fixed-point value and the dependence of i* on the initial value of /i into ( p.4| 



we find to lowest order in e that n indeed has the form (1^) with the universal number C 
given by 



C 



1 



47re 



In d = 2, M approaches logarithmically slowly. For large £ we have 



nil) 



mi 



Inserting this into (|2.4| ) we find [ plq] 
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log 
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Note the logarithmic violation of the perfect scaling of (|L6] 



(2.6) 



(2.7) 



(2.8) 



B. r>o 



We will restrict our analysis to the center of the quantum critical region in d = 2 (Fig ^ 
and d < 2 (Fig |^): the initial value of fi will therefore be fixed at /i(£ = 0) = and the initial 
value of the temperature T will be close to 0. We will only need the finite T renormalization 
group equations for /x and T which are |]17| 
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To leading order in e (leading-logs) it is sufficient to assume that for d < 2 {d = 2) u is given 
by Eqn ( p.5|) (Eqn (|2.7|) ). We will now integrate the renormalization group equations until 
a scale i* where 
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The correlation length ^ is then given by 
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while the boson density is 
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log(l - e^(^*)/^(^*)) 



(2.12) 



In the last step we have anticipated that to leading order in e it is sufficient to evaluate the 
integral directly in d = 2. 



Let us now examine the results of integrating (2^) for d < 2. To leading order in e we 



find 



^(i) = -4ee 



21 
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\2mT{e') 



Using T{i) = Te^^, it is straightforward to perform the integration and obtain from ( [^.10 ) 
the leading result for i*: 



-2r '2mkBT / I 

e = . n . „ 2e log — . 



n^A2 '-'°\2e 



From (|2.11|) we therefore deduce 
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{2e\og{l/2e)y/^ V^mksT 
and from (|2.12 ) we obtain for the boson density 
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(2.15) 



1 
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"=i^^j 4;^^°H2elog(l/2e) 



(2.16) 



These last two results are consistent with the zero-scale-factor universality of ( |1.19| ) and 
( |1.23| ) and yield properties of the scaling functions i3+_ and A^ at /i = 0. 

Finally consider properties in the quantum-critical region in li = 2. The analog of Eqn 
(ill)is 
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Integrating this to leading-log accuracy and using (|2.1CI|) we find 



„„ 2mkBT , , / h'^A'^ 
e = 4^75-TT — log log 
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(2.18) 



We therefore have from ( |2.11| ) for the correlation length 
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h 

\^2mkBT 



Uoglogih'^ Ay {2mkBT)) 
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(2.19) 



which violates the universahty of |1.19| at /i = by the double logarithms. From (p.l2|) we 
get for the boson density 
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2mkBT 1 



h' 47r 
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(2.20) 



again logarithmically violating ( 1.23| ) at /i = 0. 



III. EXACT RESULTS IN ONE DIMENSION 



We have so far determined that for small e = 2 — d the /i = T = critical field theory has 
a contact interaction of strength u* = 0{e), and all other two/multi-particle interactions 
can be neglected. Remarkably, following Haldane [jl9|, it also possible to determine the 
exact critical field theory for e = 1 or d = 1. The critical field theory then has u*\e=i = C)0 
(the bosons are thus impenetrable). Moreover, all other boson interactions continue to be 
irrelevant. In Appendix |B| we consider a one- dimensional Bose gas in the vicinity of this 
strong-coupling fixed point and demonstrate this explicitly. 

The methods of Appendix |B| and earlier works |T9| , pO| ,^ use the well-known equivalence 
between the d = 1 impenetrable Bose gas and free fermions. The field theory of the critical 
end-point at /x = T = is therefore given by the free fermion Hamiltonian 



H, 



I 



cix^^(a;) 



2m dx'^ 



fi ^f{x) 



(3.1) 



where "^p are canonical fermion fields. Correlators of this theory can only depend upon 
m, /i, and T, and the zero scale-factor universality is therefore manifest. The scaling limit 
of the correlators of the uniform magnetization can now be obtained almost trivially: the 
uniform magnetization density just measures the number of particles and its correlators are 
therefore the same as those of q^bSz'^f^f- In particular we have for the scaling function 



for the uniform magnetization in (1.23) 



M{t) 



1 



n Jo 



dy 



ey 



2-t 



+ 1 



(3.2) 



where t = ^/{kBT). This scaling function has the limiting value ^A = e*/(2-\/7r) as t -^ — oo 
in the Haldane gap region, and ^A = Vt/n as t — > oo in the Luttinger liquid region (Fig ^. 
This last result combines with ( |1.23| ) to yield ( |1.6[ ) with C = 1/tt ^. 

Observables associated with correlations of the staggered magnetization, like G^, are 
much more difficult to obtain - it is necessary to express the impenetrable Bose fields in 
terms of the Fermi fields by a continuum Jordan- Wigner transformation |^ and then eval- 
uate the correlator - a naive Wick's theorem expansion of this correlator will yield an infinite 



number of terms. Recently, Korepin and Slavnov [21 1, following earlier work of Lenard [22 
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have succeeded in resuinining this expansion and showing that all space-time dependent, 
finite temperature, correlators of the impenetrable Bose gas can be expressed in terms of 
the solution and Fredholm determinant of a linear Fredholm integral equation. Thus de- 
termination of the universal scaling function A in d = 1 has been reduced to the problem 
of solving completely an integral equation, and taking the Fourier transform of the result. 
Analytic methods can take us no further, and it is necessary to resort to numerical anal- 
ysis of the integral equations. We have begun such a numerical program, and have so far 
obtained essentially exact results for the equal-time correlations - these are described below 
in Section |111A| . It should be possible to extend our results to obtain local, time-dependent 
correlations (and hence the scaling function F in ( |1 .12] )) but we have not yet done so. A 



general picture of the form of F can be obtained from the asymptotic limits quoted in Sec- 
tion 1^; should it become experimentally useful to obtain more precise numerical results for 



F, we shall be happy to provide them. 

We also note that, recently, Korepin and collaborators [^,0] have succeeding in deter- 
mining exact results for certain asymptotic properties of G^ by applying the quantum- inverse 
scattering method to the integral equations noted above. For the equal-time G^ they ob- 
tained results for the leading and next-to-leading terms as x — > oo, while for unequal time 
correlators, both x and r were sent to cxd. Unfortunately, these asymptotic results are not 
very useful in determining experimental observables which require Fourier transformation to 
functions of momenta and frequency (the large x behavior of a function implies little about 
the small k limit of its Fourier transform). Simply Fourier transforming the asymptotic 
terms leads to results which compare very poorly with the exact results which we obtained 
by the alternative means described below. We comment on some features of these exact 
asymptotic results in Appendix 0. 

In Section [illB| we will consider the limit fi ^ ksT where it is possible to make much 



greater analytic progress in determining the scaling functions. As we have already noted, 
the lower frequency properties in this region are described by Luttinger liquid criticality, 
and it possible to use conformal invariance arguments to obtain closed-from results. 

A. Equal-time Structure Factor 

The most convenient procedure for determining equal-time correlations begins with 



Lenard's [22| result for the density matrix of the impenetrable Bose gas, which is tanta- 



mount to a formal solution of the integral equation of Ref [ETI] . His result can be written 



as [H 



G{x,T = 0-) = {0\Gf{1 - 2GF)~^\x)det{l - 2Gf) (3.3) 

where the operator Gp acts on the real axis between and x, and has the matrix elements 

/cxD rllr pik,\X X ) 

i. e. the fermion Green's function. The form (|3.3|) is amenable to rapid numerical evaluation. 
We discretize the real line between and x into A^ points, whence the operators in ( p.3|) 
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become N x N Toeplitz matrices. A straightforward extension of Levinson's algorithm ||2^ 
was then used to compute the determinant and inverse of 1 — 2Gf. The computer time 
required for this step scales only as A^^, and we were able to easily uses values upto A^ = 
10000. The results for G for large x were compared with the exact asymptotic results of 
Ref 1^, with excellent agreement. Finally, we performed a spatial Fourier transform, and 

obtained results for the scaling function B |_ of the structure factor, 5* |_(fc) defined below 

(^) and by ([T3 



Our results for some representative values of t = ^/{ksT) are shown in Fig^ A computer 

program to obtain numerical values of E ^(r, t) (r = hk/y/2'mkBT) for arbitrary r, t can be 

obtained from the authors; the accuracy is limited only by computer time, but it is possible 
to obtain 3 significant figure accuracy quite rapidly. 

We note that the zero temperature limit of S \-{k) was computed by Vaidya and 

Tracy [^]: they also pointed out that the T = 0, S |_(fc) has non-analyticities |]l9l (which 

are, however, unobservably weak) at integer multiples of 2/cf, where kp is the Fermi wavevec- 
tor of Hp 



B. Luttinger Liquid 

This is the regime where fi ^ ksT. In the fermionic description of the problem, this 
means the Fermi sea is much deeper than the temperature. In this case the ground state and 
its pertinent excitations can be described by a theory in which the spectrum is linearized 
near the Fermi points k = i/cj? determined by 

By the same token, as long as the probe frequency u and momentum k do not cause excita- 
tions that probe the sea deeply (i.e. stay in the linearized region) the Bose gas will exhibit 
characteristics of the line of finite-/i, z = 1 Luttinger-liquid critical-points at T = - the 
Luttinger-liquid scaling function for the Green's function will be denoted by Al- In terms 
of Fig ^ it means that the Al will describe the lower frequency crossover around hu ~ /c_bT. 
The ratio huj/ksT can take arbitrary values as long as both hu and fc^T remain significantly 
smaller than fi. The depth of the Fermi sea will not enter any of the calculations of Al and 
fi will enter only via the Fermi velocity c. For example, in the low density limit, c is given 
by 

-^=f^V". (3.6) 

We shall see all this happen as we analyze the exact results momentarily. The crossover 
around /la; ~ /i in Fig ^ is not part of the Luttinger hquid criticality, and is instead associated 
with the z = 2 critical end-point, and the scaling function A. It should be apparent from 
this discussion that the limits hu /ksT -^ oo and ii/ksT — ^ oo of ^ do not commute. 

Note that Luttinger-liquid criticality, and associated scaling forms, hold even when 
the condition that /i be small in (|1.14|) is violated. Suppose, however, that fc^T, Uuj and 
h'^k'^ /{2m) are much smaller than /i, and that /x itself is small so that Eqn.(^TT^) is satisfied. 
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Then zero scale-factor universality of the z = 2 critical point at fi = 0, T = must be 
simultaneously satisfied. This will lead to a compatibility condition condition between the 
scaling functions A and Al, which we shall shortly examine. 

Let us begin by writing down the scaling forms of the Luttinger liquid at T = 0. We 
know that at equal times 

G{x, r = 0") = — ; x ^ oo, T = 0, fi > 0. (3.7) 

The constant D, and the exponent t], will, in general, have a non-universal dependence 
upon the microscopic couplings. However, knowledge of D, r], and a zero-sound velocity, 
c, will universally determine all remaining hydrodynamic properties in the Luttinger liquid 
regime. For example the unequal time correlation function will have exactly the same form 
as above with x replaced by the euclidean distance va?"+~c?r^. As for r], it has a value that 
depends on the Luttinger coupling. At the point // = 0, the bosons are impenetrable and 
equivalent to free fermions. For /i > 0, the deviation from the impenetrability condition can 
be translated into a residual interaction between fermions by integrating out doubly occupied 
states (see Appendix |B]). This is the marginal coupling of the Luttinger liquid. Let us note 
for future reference that rj = 1/2 for zero Luttinger coupling. Some readers may have trouble 
reconciling this with the fact that fermion-fermion correlation functions fall as 1/x in free- 
field theory. However, we have already noted that there is a rather complicated between the 
Fermi and Bose fields and we remind the reader of the chain of transformations relating the 
two. First the hard-core bosons are described by the Pauli matrices a±. The latter are than 
converted to a single component fermions "^p by a Jordan- Wigner transformation. These 
fermions are filled upto some Fermi momentum determined by the chemical potential. When 
linearized near the Fermi points, the spinless Fermi field turns into a pair of relativistic fields 
'^L,R- The hard-core boson correlation function at equal times is 

< ^{0)^^{x) >=< ^^(0)e'"r*F(-')-i'F(-')dx'^^(a;) > . (3.8) 

If we now write 

^^(x) = ^ie-*'^^^ + ^rc''"'''' (3.9) 

and drop terms that oscillate at kf, we obtain 

< ^(O)^t(^) >^< ^^(o)e^-r(*I(-')*^(-')+*fl(^')*«(-'))'^--'^t. (3,) + L^R>. (3.10) 

Clearly this is a complicated object in the Fermi theory. To evaluate it one uses bosonization. 
Using the standard dictionary it is possible to show that it is proportional to the two-point 
function 

^ ^-iv^m^iv^k^^) >oc J- (3.11) 



where is the field dual to the usual boson field. (See for example Ref ^1 



We now consider the correlations at finite temperatures. In general the passage from zero 
to nonzero temperatures is nontrivial since in the latter case non only the ground state but 
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excited state correlations enter. However in two euclidean dimensions, in a relativistically 
invariant theory such as this one, we have the remarkable result from conformal field theory 
that was first pointed out by Cardy |2^ . Using the conformal mapping between the infinite- 
plane and the strip of finite-length Lj. = hc/ksT along the imaginary-time direction, one 
can obtain from eqn.( |3.7| ) 



Gl{x,t) 



D (2-nkBT" 



2^/2 



he 



cosh 



he 



— cos 



h 



-r,/2 



(3.12) 



(The operator e~^^^^^^ is a primary field, which allows us use conformal invariance as above.) 
The subscript L has been placed to emphasize that this formula is valid only near the lower 
frequency crossover in Fig. |^ 

The ease with we obtain this result should not detract us from noting its importance- 
rarely does one have the thermally averaged correlation functions of an interacting system. 
We shall therefore spend some time analyzing this result. 

As a first step let us extract form this the correlation function at equal time, for long 
distances. It is readily seen that 



G{x,0- 



j^f^nkBTy^^ f-r^TikBTx 



he 



he 



T > 0, fi > 0, X— >-cxD 



(3.13) 



Now let us consider the regime where (|1.14| ), is also satisfied and so zero-scale- factor 
universality holds. Thus D, rj, and c can no longer be nonuniversal, but must be universal 
functions of fi and m. Let us now invoke the asymptotic (long distance) hard-core boson 



scaling functions at the fi = critical point which are known by exact solution p6| , |23 

^/2mkBT Poo ( ^/2mkBT 



G{x,0- 



nx 



exp 



h \^t. 



(3.14) 



where t = n/{kBT). The constant poo is a known universal number; more details on this 
formula are relegated to Appendix |C|. We see that this agrees with the Luttinger liquid 



result Eqn. (|3.13|) if 



V 



1 
2 



2/i 
m 



1/2 



71 \ h^ J 



1/4 



(3.15) 



Notice also how the chemical potential entered only via the Fermi velocity as anticipated 
(excluding the /i dependence of the pre- factor D). 

We now consider the Fourier transform of ( p.l2[ ) to obtain the corresponding GL{k,iuJn) 
at the Matsubara frequencies along the imaginary frequency axis 



GL{k,iUn) 



dx dre-*('=^-'""")Gi(x,r) 

-oo Jo 



(3.16) 



Given the scaling form of Gl{x, t) it follows that the two integrals involved in the transform 
lead us to the scaling form 
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Gf{kM 



he 



2-r, 



Al 



Tiuj Tick 
ksT' ksT 



(3.17) 



where Al is a completely universal scaling function, dependent only upon the value of rj. 
We will determine Al in closed form below. 

Now if we take the limit fi/kBT — i> oo at fixed uj/kBT and hck/kBT (while satisfying 
( |1.14| ) of course), the system is described simultaneously by the Luttinger liquid result (p.l7|) 
and the zero scale-factor universality of (|1.10|) . Comparing these two results and using (p.l5|) 
we obtain immediately the compatibihty condition between the reduced scaling function Al 
at ?7 = 1/2 and the scaling function A 



AL{uJ,k)\r^=l/2 



TT 



V2po 



lim —;=A 



k 



UJ. 



2v/t' 



t 



(3.18) 



We are using here, as before, the notation ZU = huj/(kBT), k = fick/{kBT) and t = ij,/(kBT). 
We have demanded here that n enters Al only through c. Notice that the t — >■ cx) limit is 
taken at fixed cJ, and as we have noted before and shall see explicitly below, the a; — > oo 
limit of Al does not agree with the cZ; — i> oo limit of A which was implicit in ( |1.13| ) 

The remainder of this subsection is devoted to obtaining explicit results for Al and 
associated scaling functions in Minkowski space. There are two ways to proceed. 

The first approach is to perform the transform for complex (Matsubara) frequencies and 



then make the substitution iur. 



u to obtain —Giuj). The transforms are tedious to 



perform but the interpretation of the results is instructive and we shall do so soon. The 
general principles we learn about analytic continuation into the complex plane are usually 
illustrated with trivial examples (i.e., noninteracting propagators) and here we have one of 
the few nontrivial and hence instructive cases. 

The second approach is peculiar to this problem and relevant because conformal invari- 
ance methods always give the correlations in coordinate and not momentum space. Thus 
one can continue the results from imaginary to real time first and then take the transform. 
That calculation may be found in Section 3.3 of Ref pS] and has its own pedagogical value. 



Returning to the first approach, we used the identity 



X-W2 



r(V2) h 



d\\^/' 



-Ig-AX 



(3.19) 



to put the cosh and sin terms in ( |3.12|) up in the exponent. The x and r integrals were 
then analytically performed, followed finally by the A integration. After using ( ^.17| ), the 
final result gave us values of the scaling function Al at the Matsubara frequencies along the 
imaginary frequency axis. We found 



Al (iio, k) 



.'-'^(1-^ 



T] \u!\+ik\ /t] \uj\ — ik^ 
4 47r / li An 



22-r, 



r 1 



7] \uj\ + ik 
4 An 



r 1 



+ 



\UJ\ 



ik 



(3.20) 



An 



This result was obtained earlier by Schulz [^ and Shankar p9| in a different context but 
not analyzed in any detail. 
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First, the above function ^l(|cJ|) specifies our knowledge at positive and negative Mat- 
subara points. Our goal is to construct the physical real frequency correlation function and 
its singularity structure from it. As it stands, Al can only be used for numerical purposes 
and not for studying analytic structure since |cJ| is neither analytic nor ant i- analytic. In 
other words we can use Al to calculate values of the putative function at the Matsubara 
points. The mod symbol tells us that the function we are seeking has the same value at any 
postive Matsubara point and its negated image. 

Now, on the real frequency axis we have a retarded correlation function, which we assume 
is well defined. Since the factor e*'^* converges in the upper half-plane (UHP) the function on 
the real axis has an analytic extension to the UHP which is free of singularities. Its values at 
the Matsubara points uJ = 27m, with n an integer are given by Al- Such a function is readily 
found: simply drop the mod symbol on uJ in the formula for Al- Let us call the function 
(with the mod symbol dropped) ^^■'^^(a7).This A^^^ip) is the unique analytic function 
(with good behaviour in the UHP) determined by our knowledge at postive Matsubara 
frequencies. Being an analytic function it has a continuation to the lower half-plane (LHP) 
which is however not guaranteed to be free of singularities or to have anything to do with 
the original problem. In particular the poles that the F functions have in the LHP are not 
germane to the physical response function. In fact this continuation to the LHP of A^^^{u) 



does not even agree with the data we have in Eqn.( 3.20 ) for negative Matsubara points: since 
AY^^^ip) 7^ AY^^^{—uj) it is not invariant under the change of sign of frequency as the given 
data are. However, there is an analytic function which will duplicate the given data in 
the LHP: it is obtained by replacing |cZ;| by —UJ in Eqn. (|3.20|) . Such a function, A 



LHP, 



UJ 



satisifes 

A^,^^iuJ)=Ar'i-uJ)- (3.21) 



This function will agree with Ac of Eqn. ( |3.20D at points with negative Matsubara frequencies 



and be free of singularities in the LHP. However its poles in the UHP have no physical 
significance. Thus the function at real frequencies is the limit of two different functions 
as we approach the real axis from above or below. The true singularities of the physical 
response function are due to the mismatch of these two functions and not due to the poles 
they have in regions where they no longer represent the physical function or a continuation of 
it. So we must consider the difference between these two functions A^^^ip) and Al^^(P) 
on the real frequency axis. It is readily shown that the two functions are conjugates of each 
other there so that the discontinuity is just twice the imaginary part. The physical response 
function's real singularity is therefore a cut and not poles. This is reasonable since a theory 
with gapless excitations on top of the ground state (and hence a cut at zero temperature) 
cannot lose its spectral weight in these regions by the inclusion of higher states in the thermal 
average at finite temperatures. (The converse is possible: a cut free region on the real axis 
at T = can close up at T 7^ 0.) 

In Fig I we plot —ImALip, k)/in as a function of ZJ and at a representative set of values 
of k. Notice that for large k there is a well defined peak at a; ~ fc: this is a signature of the 
propagating modes in the Luttinger liquid ground state and represents the behavior of the 
intermediate region ksT ^ cu ^ /i in Fig ^ At smaller values of k notice that the peak in 
Fig ^ remains at oJ = 0. This is the z = 1 quantum relaxational behavior (Fig ^ where the 



strong interaction between the thermally excited Luttinger modes has left only overdamped 
excitations. 

We now consider a couple of other experimental observables, related to local and equal- 
time correlations respectively. 

1. Local Green's function 



The local Green's function, Gf was defined in ( |1.11|) . In the Luttinger liquid regime, we 
can deduce that, provided cu -C /i, this observable satisfies the scaling form 



ImGf^icu) 



-sgn cj — \ujr Ft - — — 



(3.22) 



where F^ is a universal function, specified completely by the function Al in ( |1.11| ). Moreover, 
as the Luttinger-liquid criticality has a particle-hole symmetry, F^ must be an even, positive 
function of cJ. As already noted, we expect on general grounds that ImG^(co') ~ uo for small 
uj at finite T: therefore Fl{uJ) ~ |^P~^ at small ZJ. We also note that for rj > 1 the real part 
of the local Green's function will not satisfy an analogous because the integral in ( p. . 1 1 ) is 
then dominated by large momentum contributions. 

There is again a compatibility condition between the Luttinger liquid scaling function 
Fl and the scaling function F in ( |1.12| ) quite analogous to that for Al, A in ( |3.18| ); we have 



Fr(oj) 



IT 



|r?=l/2 



lim FCuJ, t) 



(3.23) 



As before, the limits U ^ oo and t ^ oo do not commute, and the To ^ oo limit of the 
exact Fl computed below will not agree with that of F in (|1.13|) 

Let us finally present the exact computation of F^. We use the result ( ^.121 ) at x = 0, 
Fourier transform to Matsubara frequencies, analytically continue and take the imaginary 
part to obtain the following result for F^ 



Fl(P 



7^\ — IttII"'' 



UJ\ 



TT 



"-1/2 gjnh 




A plot of this function is shown in Fig | 
For small aJ we have 



FLiuJ) 



2r 



Y !l_ ^1^1 



1 + r] 



27r 



\UJ\ 



r(^ 



(3.24) 



(3.25) 



This is the behavior characteristic of the z = 1 quantum-relaxational regime of Fig ^ It is 
expected that the limits uJ —>■ and t ^ oo do commute, so combined with ( |3.23| ), the above 
result gives us the small uJ behavior of F{uj, t) at large values of t. 
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In the opposite limit of large uo we crossover to the critical correlations of the Luttinger 
liquid ground state in which case 



Fl{^) 



2i-%3/2 



r(l±^k^^ 



\UJ\ 



oo 



(3.26) 



This last result can also be obtained by a Fourier transform of the relativistic zero temper- 
ature correlator. 



2. Structure Factor 



The two structure factors S^-{k) and 5.+ of the antiferromagnet were defined in Eqn 
( |1.18| ). The Luttinger liquid behavior has particle-hole symmetry so in this regime, the two 
structure factors are essentially equal and will be denoted by by the common value Silk). 
The scaling form for 5'l(A;) follows from ( |1.16| ), ( [L.18|) and the scaling of Gf in (|3.17|) : 



S^k) = ZD 



he 



i-v 



Br 



Tick ' 
k^, 



(3.27) 



where the constant Z was introduced in ( |1.16| ) and Bl is a universal function obtained below. 
There is a compatibility condition between the Luttinger liquid scaling function Bi and the 
scaling function i3_|__ in ( |1.19D which is quite analogous to that for Al-, -A in ( |3.18| ): 



BMU-. 



TT 



1/2 



2po 



lim -^B<^ 



k 



*^°° y/i 



2Vt' 



(3.28) 



Using the other scaling function B ^ on the right-hand side would yield, from (|L2CI|) , an 

identical result. Again the limits t ^ oo and k ^ oo are not expected to commute. 

Finally, the exact determination of Bl'. we simply perform a spatial Fourier transform 
of ( p.l2|) at r = 0; we obtain in this manner 



BL{k) 



T 



TT 



■n-i/2_ 



Hhi 



r i + lt 

V2 2n, 



(3.29) 



This result is illustrated in Fig |^. For small k, Bl reaches a constant whose value is easily 
obtainable from (|3.29| ). For large k we find from (|3.29|) 



2i-vi/2r 



BL{k) 



1 — rj 



t-" ' 



\k\ 



oo 



(3.30) 



Again this last result could have also been obtained by direct computation at T = 0. 
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IV. CONCLUSIONS 

This paper has studied the universal, finite temperature properties of a dilute Bose gas 
with repulsive interactions in dimensions less than or equal to 2. In the vicinity of the T = 
onset at zero chemical potential, /i, it was argued that the leading scaling properties obey, 
for d < 2, a. hypothesis of zero scale- factor universality. This means that the entire two-point 
correlator is a universal function of just /i, T and the bare boson mass m. 

The main motivation behind this study is the mapping onto it of the properties of 
quantum-disordered antiferromagnets in a finite field. In particular, in d = 1, Haldane gap 
antiferromagnets undergo a magnetization onset at a critical field which is expected to be 
in the universality class of the Bose gas transition. Applicability of our theory requires that 
there be no spin anisotropy in the plane perpendicular to the applied field. Most materials 
do have some anisotropy - in this case we would require that the temperature, T, be larger 
than the anisotropy gap, before applying our results. 

We have described the rather complicated properties of numerous scaling functions, which 
may rather difficult to disentangle experimentally. A useful starting point for neutron scat- 
tering experiments appears to be the following. Perform the experiment somewhere in the 
Luttinger liquid region where the absolute value of the scattering cross-section is also the 
largest. Measure the local susceptibility Gf{u) and see if collapses onto the scaling form 
( |1.12| ). For large fi/ksT, we have a rather complete picture of the scaling function F: for u 
smaller than or around kgT, we can deduce F from (|3.23|) and (|3.24|) , while for extremely 
large u we can use ( |1.13| ). 

Another possible application of our results may be to quantum-disordered antiferromag- 
nets in d = 2. By measuring the ground-state magnetization in a field, and comparing the 
result to (pT^), ( p..21| ) and ( p.8|) it may be possible to determine the spin S^ of the elemen- 
tary excitations above the ground state. Of course, we would also need an independent 
determination of the quasi-particle mass m. 
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APPENDIX A: MAGNON OPERATORS 

It is clear that the predominant coupling of neutrons will be to the antiferromagnetic 
order parameter (p^, with a = x,y,z. In a zero-field spin-fiuid phase with confined spinous, 
this order parameter corresponds to a real, massive, bosonic triplet. In this appendix, we 
want to explore in some more detail the relationship between (p^ and the complex bosonic 
field, \l/ in in Eqn. (|1.3 | ) . We note that the field theory (|1.3|) will also describe the magneti- 
zation onset transition in antiferromagnets with deconfined spinous; however, in this case, 
the relationship between the neutron scattering cross-section and the field \1/ will be quite 
different |T^, and will not be considered in this paper explicitly. 
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For simplicity, we consider d = 1, although the analysis is quite general. First we expand 
the real triplet, 0^, in terms of magnon creation and destruction operators as usual: 

/oo fjU 1 
_^^^/^Mk)e''^ + aie-^'^] (Al) 

where 



Uk = VA2 + P, (A2) 

A is the Haldane gap, and a^, a)^ ^^^ magnon destruction and creation operators for each 
of three polarizations. Let us focus on the two combinations 

0±(x) = ^^ = /_^ ^-^[(a.(A:) ±m,(A;))e^^- + {al{k) ± tal{k))e-^'% (A3) 

Observe that (f)± are adjoints of each other, but commute with each other. Let us now 
consider an effective theory for energies far below the Haldane gap. In this case we can 
make the replacement 



and obtain 

0±(x) 



cc'fc ^ A (A4) 



flk ^ 

-.[{a,{k) ± mj,(A;))e*'^" + {a\{k) t ia\{k))e-''^'^] (A5) 



-oo 27r V4A 
^ [^>±{x) + ¥^{x)] (A6) 



/2A 



where \l/±(x) destroys a spin ihl magnon at x while \l/^(x) creates a spin =i=l magnon at x. 
Suppose we next argue that, when the applied uniform field is near its critical value, only 
the spin up magnon (very light) will either be easily created or destroyed, so that we may 
drop the spin down creation and destruction operator in the above expressions. Then we 
obtain 

^ =*+(^) (A7) 



/2A 

-.'^\{x). (A8) 



V2A 

Observe that now (up to a scale factor) 0-|-, which were previously commuting, are now 
canonically conjugate fields. This is just like in the Hall effect wherein x and y , which are 
commuting coordinates in the full Hilbert space become conjugates in the lowest Landau 
level. It is also clear from the discussion that the field \1/ in the coherent space integral in 
Eqn. ([T7^ ) is precisely this complex conjugate pair. 
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APPENDIX B: BOSON HUBBARD MODEL IN ONE DIMENSION 

Consider bosons 6, moving on the sites, i of a chain described by the Hamiltonian 

H = -wY. {blh+i + bl^.h - 2blh) + Y. {Vn,{n, - 1) - ^n,) (Bl) 



where rii = h\hi is the number operator, w is the hopping matrix element, and V is the 
on-site repulsion between the bosons. In the limit of large V , states with more than one 
boson on a site will only occur rarely, and it should pay to restrict the Hilbert space by 
projecting out such states. However, the elimination will generate a residual interaction of 
order w"^ /V between the states on the restricted space. This interaction can be determined 
by the usual second-order perturbation theory and leads to the effective Hamiltonian 



H. 



-w 



i i 

2w^ 
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Y, {2blbl^,bi+ibi + blbUh+ibi + blbl^.k.ibi) 



(B2) 



We reiterate that Hf, is non-zero only on states with at most one boson per site. Notice now 
that this reduced Hilbert space is identical to that of spinless fermions. The transformation 
between the bi and the spinless fermion operators fi is of course the Jordan- Wigner mapping 



k = n(i - 2/1/,)/, 



(B3) 



j<i 



We now insert ( P3|) in ( [B2|) and take the continuum limit with fi = ^/a'^F{x 
h /(2?Tia^) (a is the lattice spacing) and obtain 



«a , w 



H, 



dx 



vl> 
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V dx 
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(B4) 



It is now clear by power counting that the four-Fermi coupling term is clearly an irrelevant 
perturbation to the /i = T = fixed point. It is in fact also not difficult to see that all inter- 
actions between the fermions are irrelevant, and that this result is not special to the Boson 
Hubbard model considered here. The key point is of course that a term like '^p'^p^p^p, 
which is the only interaction term which is relevant by power counting about the free fermion 
fixed point at /x = 0, vanishes identically because of the fermion anticommutation relations. 
The significance of the four- Fermi coupling changes when we consider the scaling dimen- 
sions of operators about the Luttinger liquid fixed points. In this case we decompose the 
fermion field into left- (^l) and right- (^_r) moving excitations with a linear dispersion, and 
obtain the long-wavelength Hamiltonian 



Hr 



dx 



he \E' 



t/^ 
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where c = Tikp/m and the Fermi wavevector kp is given by h'^kp/(2Tn) = /i. Performing 
power-couting on the z = 1 free field part of Hp we now find that the four- Fermi coupling is 
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now marginal. Note however that the coefficient of this four- Fermi couphng is suppressed by 
a factor of kp, which vanishes as one approaches the z = 2 critical end-point. By the usual 
logarithmic perturbation theory at T = we can determine that the four-fermi interaction 
modifies the exponent rj of Section [1IIB| by 



V = ^-T 



'W/J, 



V 



(B6) 



where 7 is a numerical constant of order unity. Notice, as expected, that the correction to 
7] vanishes as /i — >■ 0. It is also apparent that the impenetrable limit ^ — >■ 00 is equivalent 
to the vanishing density (/i — > 0) limit. 

At finite T, there will be corrections to the correlators with terms like 
{^/tfI/V)log{f^/kBT). These cannot be neglected when 



ksT < n exp 



V 



-7 



/w/i 



(B7) 



where 7' is of order unity. Zero scale-factor universality is thus violated for arbitrarily small 
H, when T is smaller still and satisfies (p7|) . Notice however that the boundary specified by 
( [B7D lies well below ksT ~ /z crossover to the Luttinger liquid regime (Fig|l|). 



APPENDIX C: IMPENETRABLE BOSE GAS IN ONE DIMENSION 

Its et. al. [^] have recently obtained some exact asymptotic results for the equal-time 
boson Green's function of the d = 1 impenetrable Bose gas. Recall that this model is 
precisely the scaling limit describing the z = 2 quantum phase transition with zero scale- 
factor universality. In this Appendix, we show how the requirement that the scaling functions 
be analytic in fi/ksT, can lead to a considerable simplification of their results. 

The asymptotic results of Its et. al. can be written in the form 



G{x, 
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exp 
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cLS X 



CXD. 



(CI) 



where A{t) and f{t) are functions to be determined (as before t = ^/ksT). From the 
arguments in Section |T|, it is clear that A{t) and f{t) are also universal crossover functions 
of the /i = 0, T = 0, quantum phase transition the repulsive, d = 1 Bose gas with arbitrary, 
short-range interactions. Its et. al. obtained two separate, closed-form, integral expressions 
for f±{t) and A±{t) valid respectively for t > and t < 0. The two expressions were quite 
distinct and there appeared to be no straightforward relationship between them. 

Here, we point out that the absence of any singularity in the finite T Bose gas, in fact 
requires that f{t) and A{t) be analytic for all finite, real values of t. In other words, the 
functions /+(t) and /_(t) must be analytic continuations of each other (similarly for A^{t) 
and A_(t)). We have in fact succeeded in proving that the expression of Its et. al. for 
f+{t) is in fact the analytic continuation of their result for f-{t). We have been unable to 
establish a similar result for A± (t) , but have performed numerical tests on their expressions, 
which leave essentially no doubt that A is also analytic. 
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With the help of the above considerations, it is possible to deduce from Ref p3| a simple 
closed- form result for f{t) and A{t) which is valid for all t 
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The value of the constant poo was obtained by matching to the T 
Tracy |2^ 



(C2) 



result of Vaidya and 
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0.924182203782.... 



(C3) 



Ag being Glaisher's constant [^. Note that the analyticity of / and A for all real t is 
manifest. We have plotted the functions f(t) and A{t) in Fig H. They obey the asymptotic 
limits 
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Both asymptotic limits of f{t), and the t — >■ +oo limit of A{t), can be obtained directly from 
(|n2|). The t —>■ — oo limit of A{t) is more difficult to obtain from (C2), and we used instead 
the second expression for A{t) in Ref |2^. Demanding that these two methods of obtaining 
the limit be identical in fact provides one with an independent derivation of the value of the 
constant poo ! 

We also recall [^ that 



/(O) 



C(3/2) 



1- 



n 



2V2, 



0.95278147061075... 



(C5) 



We have already noted that taking a Fourier transform of the asymptotic results (|C1D 
to obtain the structure factor, yields results which compare very poorly with numerically 
exact results of Fig ^. 

Finally, we note that the requirement of analyticity as a function of t should apply to 



essentially all of the equal-time and unequal-time results of Its et. al. [p3| , p4| . In every case, 
they have obtained separate expressions for t < and t > 0: proving that these are analytic 
continuations of each other will lead to highly non-trivial checks on the results, and should 
also produce some fascinating mathematical identities. 
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FIGURES 
FIG. 1. Phase diagram in d = 1. The dashed hnes indicate crossovers. 

FIG. 2. Phase diagram in d = 2. The dashed line is a crossover while the full line is a 
Kosterlitz-Thouless phase transition. The field H has absorbed a factor of g^B- The location of 
the Kosterlitz-Thouless transition is determined from Refs |1(:,17| 



FIG. 3. Properties of the different regimes of Fig |l| and |2| as a function of the measurement 
frequency uj (the wavevector A: ~ 0). All crossovers are described by the universal scaling function 
A. The crossover for /i ^ ksT in d = 1 near huj ~ ksT is also described by the Luttinger liquid 
scaling function Al 

FIG. 4. Exact results for the scaling function B \-{r,t) for the structure factor (Eqn ( |1.19| )) in 

d = 1. We have r = hk/\/2mkBT, and t = fj,/{kBT). We have chosen some representative some 

representative values of t; a computer program to evaluate B |- for arbitrary t is available from 

the authors 

FIG. 5. Exact values of the scaling function — Im^i(U, A;)/aJ (given in Eqn ( ^.20D ) for the 



Green's function in d = 1 in the Luttinger liquid regime at ry = 1/2 - Eqn (3.17). We have 



uJ = huj/{kBT), k = hck/{kBT). The values for k are 3 times larger than those on the graph. 
Notice how the spectrum evolves from an overdamped, relaxational peak at small k (as for k = 1) 
to a damped, propagating mode at large k (as for A; = 3, 5). 



FIG. 6. Exact scaling function Fl{uj) (given in Eqn ( 3.24[ )) for the imaginary part of the 



local susceptibility in d = 1 in the Luttinger liquid regime at 77 = 1/2 - Eqn ( 3.22| ). We have 
U = huj/{kBT). 



FIG. 7. Exact scaling function Biik) (given in Eqn ( 3.29| )) for the structure factor in d = 1 in 



the Luttinger liquid regime at ry = 1/2 - Eqn ( 3.27| ). We have k = hck/{kBT). 



FIG. 8. The scaling functions f{t) and A{t) {t = ii/{kBT)) of the d = 1 Bose gas defined by 
Eqns dClD and (|C^) 
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